Metal-insulator transitions: 
Influence of lattice structure, Jahn- Teller effect, and Hund's rule coupling 
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We study the influence of the lattice structure, the Jahn- Teller efi'ect and the Hund's rule coupling 
on a metal-insulator transition in AnCeo (A= K, Rb). The difference in lattice structure favors 
A3 Ceo (fee) being a metal and A4C60 (bet) being an insulator, and the coupling to Hg Jahn- Teller 
phonons favors A4C60 being nonmagnetic. The coupling to Hg (Ag) phonons decreases (increases) 
the value Uc of the Coulomb integral at which the metal-insulator transition occurs. There is an 
important partial cancellation between the Jahn- Teller effect and the Hund's rule coupling. 



The competition between the Coulomb repulsion, the 
kinetic energy, the Jahn- Teller effect and the Hund's 
rule coupling leads to interesting physics. Examples are 
perovskites, e.g., the manganites [||, and alkali-doped 
fuUerenes [Ej . Here we focus on the metal-insulator tran- 
sition for an integer number of electrons per site. This is 
particularly relevant for the fullerenes, since A3C60 (A= 
K, Rb) is a metal ||] while A4C60 is a nonmagnetic in- 
sulator ^Q . According to band theory both are metals 
1^, and A4C60 must therefore be an insulator due to in- 
teractions left out in band structure calculations. 

The metal insulator transition in a correlated system is 
usually discussed in terms of the ratio U/W R], where U 
is the Coulomb interaction between two electrons on the 
same molecule and W is the one-particle band width W. 
The ratio U/W is, however, almost identical for AaCgo 
and A4C60 |§|j||J|]- The question is then why not both 
systems are either metals or insulators. To study this, 
we apply the dynamical mean- field theory (DMFT), pro- 
jection Quantum Monte-Carlo (QMC) and exact diago- 
nalization techniques to models of A„C6o- 

For the Fullerenes it is believed that U/W ^^ 1.5 — 2.5 
. In spite of this large ratio, these systems are close 
to a metal- insulator transition due to the orbital degen- 
eracy TV = 3 of the partly occupied ii„ band |||,|l2l ■ The 
lattice structure is fee for A3C60 and bet for A4C6o- The 
important electron-phonon coupling is to Hg Jahn- Teller 
phonons. We find that the difference in lattice structure 
alone can explain why A3C60 is a metal but A4C60 is 
an insulator and that the electron-phonon coupling can 
explain why A4C60 is nonmagnetic. We find an impor- 
tant competition between the Jahn- Teller effect and the 
Hund's rule coupling. The Hg and Ag intramolecular 
phonons are found to have the opposite effect on the crit- 
ical Uc, ioT which the metal-insulator transition occurs. 

We consider a model of A„C6o which includes a three- 
fold degenerate ii„ level on each molecule and the hop- 
ping between different molecules 



where ijlam creates an electron on molecule i with the 



quantum number m and spin a. The hopping matrix el- 
ements tijmm' pil include the orientational disorder |12| ] 
and the lattice structure, with nearest neighbor hopping 
for the fee structure and a weak second nearest neighbor 
hopping for the bet structure Q. The Coulomb interac- 
tion is given by 



H\S — Uxx 2_^ nira^Uiml + Uxy 2_^ / , 



iaa m<m 



+ 2^^^ ^ ^irnAlm'a'^^na'Arn'a (2) 

icrcr m^m 

where U^x and Uxy describe the interaction between 
equal and unequal orbitals, respectively. K is an ex- 
change integral and Uxx — Uxy + 2K. Finally we include 
the interaction with a five-fold degenerate Hg phonon on 
each site 



ifph = i^ph J2 ^l 



E 






(brMu), 



where bi^, creates a phonon with the quantum number 
v on the molecule i. The matrices V , are deter- 

A — h ■ rnrn 

mined by symmetry |13| . The coupling constant g is 
related to the dimensionless electron-phonon coupling 
A ~ {5/3)N{0)g'^/ujph- We also consider the couphng 
to A„ phonons, for which V^ , is diagonal in m and m . 

y ^ ' mm ^ 

In a first step we analyze the effect of the lattice 
structure alone, neglecting the electron-phonon coupling 
{g = 0) and the multiplet effects {K = and Uxx = 



Ux 



U). We use a projection Quantum Monte-Carlo 



(QMC) r = method in the fixed node approximation 
Pjl, which gives quite accurate ground-state results for 
this model P]. A3 Ceo and A4C60 differ in the number n 
of conduction electrons per site and in the lattice struc- 
tures. For a fee lattice, n = 3 and n = 4 give Mott 
transitions at almost the same U^ llql. We therefore fo- 
cus on the difference in lattice structure, and consider 
n = 4 for clusters with M molecules put on fee or bet 
lattices. The band gap for filling n is 



Eg = E{nM + 1) + E{nM - 1) - 2E{nM), 



(3) 



where E{N) is the energy of a system with N electrons. 
We want to extrapolate to M — > oo and determine the 
Uc for which Eg is zero. To reduce the finite size effects 



we add 



Eg{U) = Eg{U) 



U 
M 



Eg{U = 0), 



(4) 



where Eg{U = 0) is the band gap for U = 0. These 
corrections go to zero for large M, but they improve the 
extrapolation A/ ^ oo. Fig. [^ shows that the metal- 
insulator transition happens for a substantially smaller 
U/W for the bet {Uc/W ~ 1.3) than for the fee structure 
{Uc/W ^ 2.3). The insulating state is antiferromagnetic. 
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FIG. 1. The energy gap Eg (Eq. (0)) as a function of 
U/W for M = 32 molecules on fee and bet lattices and the 
filling n — 4. The lines are guides for the eye. 

To understand these results, we note that on the fee 
lattice it is possible to hop on a triangle, i.e., to return 
to the original site after three hops. On a bet lattice, on 
the other hand, this is not possible if the small second 
nearest neighbor hopping integrals are neglected. The 
simplest systems with these properties are a triangle and 
a square, each site having a level with spin but no orbital 
degeneracy. A nearest neighbor hopping integral t < 
connects the orbitals. The one-particle spectrum is ±2i 
for the square and —2\t\ and t for the triangle. For the 
triangle there is a state with maximum bonding character 
(— 2|t|), but it is not possible to construct an optimally 
anti-bonding state, due to the presence of frustration. 
Thus the one-particle band width are W — 3t and 4t for 
the triangle and the square, respectively. The curves in 
Fig. P mainly differ in the large U limit and we there- 
fore consider this limit. We construct the many-body 
states of the triangle with two, three and four electrons, 
which determine the band gap (Eq. dq)). The energy 
E{3) — 0{t'^/U), since hopping is suppressed to order 
t/U. For the case of four electrons, we construct all states 
with the minimum (one) double occupancy and S^ = 0- 
These states describe how the double occupancy hops 
around the triangle. The original state is, however, not 
recovered after one loop, since the spins on the sites with 
a single occupancy have been flipped. Moving the double 



occupancy around the triangle a second time restores the 
spins and the original state is recovered after six moves. 
The corresponding 6x6 matrix has the extreme eigenval- 
ues ±2i. In the lowest many-body state of the triangle 
with four electrons, it is therefore not possible to restore 
the state in an odd number of hops, and the frustration 
does not show up. In a similar way we obtain the lowest 
energy — 2i for the two-electron state. The square has 
the same energies. Thus 



Eg = U-'i\t\ = U-^W 
Eg = U -A\t\ = U -W 



for a triangle 
for a square 



(5) 



Both the triangle and the square have no frustration in 
their many-body states, and for fixed t the gaps are the 
same. The one-particle band width W, however, is re- 
duced by the frustration in the triangle, and expressing 
the Eg in terms of W requires a larger prefactor in the 
frustrated case. These results give a qualitative explana- 
tion of Fig. ||. 

Although the calculation above can explain why A4C60 
is an insulator while A3C60 a metal, it incorrectly pre- 
dicts A4C60 to be antiferromagnetic. The calculation ne- 
glects, however, the coupling to the Jahn- Teller phonons, 
which tends to make A4C60 a nonmagnetic insulator p6[ . 
The electron-phonon interaction has been estimated from 
photoemission experiments for a free molecule |1^. We 
describe the eight Hg phonons by an effective mode, 
with the logarithmically averaged frequency Uph = 0.089 
eV, and the effective coupling g = 0.089 eV. For a free 
molecule this leads to a singlet being 0.29 eV below the 
lowest triplet. This triplet-singlet splitting is larger than 
an experimental estimate of 0.1 cV for A4C60 0- The 
splitting is, however, reduced by the competition with 
the Hund's rule coupling. An estimate of the exchange 
integral K based on an ab initio SCF calculation |lq ] 
gave K = 0.11 eV ||l^. This number is, however, ex- 
pected to be reduced by correlation effects. For instance, 
for atomic multiplets a reduction by 25 % has been found 
P0| . Indeed, we find that the experimental triplet-singlet 
splitting is reproduced by using K = 0.07 eV. 

Since the metal-insulator transition depends on a com- 
petition between the kinetic and Coulomb energies, and 
since we may expect the electron-phonon coupling to 
reduce the hopping, we may expect this to reduce Uc 
pl| . We therefore study the effect of phonons on Uc (for 
K = 0). 

For this purpose we apply the dynamical-mean field 
theory DMFT [Q. We use hopping integrals for a 
Bethe lattice in the infinite dimensional limit t^^.^i ~ 
t*S^^i /^/z, where z — > cx3 is the connectivity. The im- 
purity model, resulting in the DMFT, is solved with a 
QMC method p2|. The phonon fields are treated fully 
quantum mechanically, and they are updated together 
with the Fermion auxiliary fields in each Monte Carlo 
step p3[. We use the one-particle band width W ~ 2 
and a Trotter break up At = 1/3. 



For an insulator G{t = /3/2) decays exponentially with 
/3, where G{t) is the electron Green's function on the 
imaginary time axis. We therefore use G(/3/2) to deter- 
mine whether the system is a metal or an insulator. 

We first compare the coupling to Ag and H^ phonons 
for 71 = 3. Fig. ||a shows that G(/3/2) is reduced as 
U/W is increased, since the system gets closer to a metal- 
insulator transition. For A = extrapolation suggests a 
rather large Uc/W. For Hg phonons an increase in A 
leads to a rapid reduction of G{/3/2) and Uc, while for 
Ag phonons this leads to an increase in G(/3/2) and Uc- 

To understand these results we study a free molecule 
(Table |) and a system consisting of two molecules 
(dimer) (Table ^) in the limit 

K ^ -^ = EjT <t: LUph < T^ < C/. (6) 

Table H shows the energy gap of the dimer. In agree- 
ment with the fuU DMFT results {K = and n = 3) the 
gap is increased by a coupling to H^ but decreased by a 
couphng to Ag phonons. We first consider the Ag case. 
Since V^„j„/ = S^m' ^rnm wc cau transform the electron- 
phonon coupling to the form 



9^i^M 



n){h + blh 



(7) 



where n.i is the total occupation number operator for site 
i and n is the (integer) filling. An irrelevant constant has 
been neglected. We first study the state with 2n elec- 
trons. In the limit W <^ U hopping is suppressed, and 
m — n « 0. The coupling (Eq. (^) is then negligible, and 
the electron-phonon contribution to the energy is small. 
In the case of an extra electron or hole, however, this ad- 
ditional charge can hop even for W <^U . The coupling 
to the phonons then lowers the energy, and according to 
Eq. (^ this reduces the gap. 

For coupling to Hg phonons, the state with 2n elec- 
trons can lower its energy via the (dynamic) Jahn- Teller 
effect. Since hopping is very efficiently suppressed, the 
energy gain is accurately given as twice the energy for 
a free molecule (Table |). In the case of an extra elec- 
tron or hole, on the other hand, hopping dominates over 
the Jahn- Teller effect in the limit (y). The system can 
then only take advantage of this effect to the extent that 
it does not interfere with the hopping. The electron- 
phonon coupling then gives a much smaller lowering of 
the energy than for the state with 2n electrons, which 
increases the gap (Eq. (Bl)). 

Fig. pp shows results for coupling to Hg phonons and 
filling n — 4. Uc/W is smaller than for n = 3, although 
the lattice structure is the same as for n = 3. This can 
be understood from Table I, which shows that the en- 
ergy gain in the free molecule due to the electron-phonon 
coupling is larger for n = A. This enters in E{nM), 
while the electron-phonon coupling plays a smaller role 



for E{nM ± 1). The electron-phonon coupling alone 
would then tend to favor A4C60 being an insulator and 
A3C60 being a metal. As we will see below, this effect is, 
however, partly cancelled by the Hund's rule coupling. 
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FIG. 2. (a) The electron Green's function G{t = /3/2) as 
a function of U/W for the filling n = 3 and different values 
of the electron-phonon coupling A. The figures compares the 
coupling to Ag and Hg phonons. (b) G(/3/2) for coupling to 
Hg phonons and n = 4. G(/3/2) ~ implies an insulator. 

The coupling to the Hg phonons pushes Uc for A3 Ceo 
to the lower end of the physical range of U/W, raising 
some questions of why not also AsCgo is an insulator. Al- 
though, the Ag phonons tends to increase Uc, this should 
not be important due to the weak coupling to the Ag 
phonons [M. However, there is a substantial coupling to 
a plasmon in A3C60 PH]- This should tend to increase 
Uc, since it couples to the electrons in the same way as 
the Ag phonons. Below we show that the Hund's rule 
coupling also plays an important role in this context. 

We next consider the effects of the Hund's rule cou- 
pling {K > 0). Since these terms in Eq. (||) lead to 
a sign-problem in the DMFT QMC calculation, we use 
exact diagonalization. To reduce the size of the Hilbert 
space we consider a four-site system with two-fold orbital 
and phonon degeneracies. The nearest neighbor hopping 
^im,jm' = ^ij^mm '^^ choseu randomly, thus reducing the 
degeneracy and the one-particle spacing. We limit the 
size of the Hilbert space by allowing a maximum of two 
phonons per site. Due to this limitation, the calculation 
is not fully converged for the larger coupling constants 
considered below. From the finite size corrected band gap 
Eg{Uxx) we estimate the critical Uxx as Uxx — Eg(Uxx), 
shown in Fig. B[ The figure illustrates that for A = an 
increase in K leads to a decrease in Uc flQl. In analogy 
to the discussion for the Jahn- Teller effect, the Hund's 
rule coupling can effectively lower the energy of the state 
with nM electrons while for the states with nAI ± I elec- 
trons, the stronger interference with hopping leads to a 
smaller lowering of the energy. For A > the competi- 



tion between the Jahn- Teller effect and the Hand's rule 
coupling tends to reduce the influence of either effect on 
Uc- This is shown in Table | and |l| and in Fig. |[ 
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FIG. 3. The estimate Uxx — Eg of the critical Uxx as a 
function of A = 4,0^ / {uiphW) for different values of K and 
Uxx/W = 3. 

To summarize, we have found that the difference in 
lattice structure favors A3C60 being a metal and A4C60 
being an insulator. The Jahn- Teller effect wins over the 
Hund's rule coupling, making A4C60 a nonmagnetic insu- 
lator. The couphng to the Hg phonons tends to strongly 
reduce the critical U for a metal-insulator transition, rais- 
ing questions about why not also A3 Ceo is an insulator. 
This effect is, however, partially cancelled by the Hund's 
rule coupling. The coupling to plasmons tends to further 
increase the critical U . 

This work has been supported by the Max-Planck- 
Forschungspreis . 
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TABLE I. The ground-state energy E{n) of an isolated 
molecule for n electrons. The quantity n(n — l)Uxy/2 has 
been subtracted. The results are symmetric around n = 3. 

n E{N) = E{N)-N{N -l)Uxy/2 



Low spin {K < ^Ej 



High spin {K > fgjT) 



]EjT 



-IQEjT+AK 



K 



-3K 



TABLE II. Eg{n) — Uxy — d-i{n)t for a two-site model as 
a function of the filling n. The hopping contribution to the 
gap is d(n)t, where d3(n) = -3, -5 and -6 for n = 1, 2 and 
3, respectively. The coupling is to one Ag (A) or one Hg (H) 
phonon per site. The results are symmetric around n = 3. 
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Eg{n) — Uxy — dz{n)t 
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